The fate of topological-insulator surface states under strong disorder 
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Three-dimensional topological insulators feature Dirac-like surface states which are topologically 
protected against the influence of weak quenched disorder. Here we investigate the effect of surface 
disorder beyond the weak-disorder limit using large-scale numerical simulations. We find two qual- 
itatively distinct regimes: Moderate disorder destroys the Dirac cone and induces diffusive metallic 
behavior at the surface. Even more remarkably, for strong surface disorder a Dirac cone reappears, 
as new weakly disordered "surface" states emerge in the sample beneath the disordered surface layer, 
which can be understood in terms of an interface between a topological and an Anderson insulator. 
Together, this demonstrates the drastic effect of disorder on topological surface states, which cannot 
be captured within effective two-dimensional models for the surface states alone. 



Topological insulators (TIs) are a novel form of insu- 
lators, where topological properties of the band struc- 
ture lead to the formation of metallic states at the sur- 
face [l|, |2|. Two-dimensional (2d) TIs were predicted 
3,1^ and then realized in quantum- well heterostructures 
5|. Subsequently, three-dimensional (3d) generalizations 
were discussed theoretically [6|-l8| and realized afterwards, 
primarily in Bi-based compounds [9|,ll0|. 3d TIs are char- 
acterized by a set of topological quantum numbers, which 
lead to a classification into strong and weak TIs. 

Strong TIs possess low-energy surface states which can 
be described in terms of massless 2d Dirac electrons, with 
an odd number of Dirac cones per surface. Remarkably, 
these surface states are topologically protected against 
the influence of disorder due to their helical nature: The 
spin-momentum locking of the Dirac theory implies that 
the states k and —k have opposite spin. This in turn sup- 
presses the lowest-order matrix element for backscatter- 
ing, k ^ —k^ leading to the absence of weak localization 
for weak potential disorder. What about strong disor- 
der? Theoretical results obtained within 2d continuum 
Dirac models suggest that such states remain metallic 
even for arbitrarily large disorder: The bet a- function de- 
scribing the renormalization-group flow of the conductiv- 
ity is always positive llMl3|. However, by construction 
the applicability of these effective theories to TI surface 
states is restricted to energies much smaller than the bulk 
gap. A, and the effect of surface disorder with a strength 
7 > A can therefore not be captured. Intuitively, this is 
because strong impurities will influence both surface and 
bulk states, causing a non-trivial coupling between the 
two which cannot be described by a theory for surface 
states alone. 

To the best of our knowledge, the effect of strong sur- 
face disorder has only been discussed locally for isolated 
impurities [1^, but the global fate of TI surface states 
under moderate or strong disorder is not known. This 
issue is not only of fundamental interest in order to un- 



derstand the range and limitations of "topological pro- 
tection" , but is also of enormous practical relevance: Im- 
purity potentials can easily be in the energy range of 
1 eV, thereby exceeding typical band gaps of 3d TIs, e.g. 
0.3 eV for Bi2Se3. Moreover, surfaces are often intrinsi- 
cally disordered, and in addition can be easily disordered 
intentionally to perform controlled studies of disorder ef- 
fects. Therefore, non-perturbative theoretical studies of 
surface disorder are clearly called for. 

In this paper, we present a comprehensive study of the 
effect of surface disorder on 3d TIs by means of large-scale 
numerical simulations of a 3d lattice model. Our central 
results can be summarized as follows: Topological pro- 
tection is only operative for disorder strengths which are 
small compared to the bulk gap, 7 <C A. For moderate 
disorder, 7 ~ A, the surface layer instead shows diffu- 
sive metallic behavior, and the momentum-space Dirac- 
like structure of its states is lost. Moreover, low-energy 
states start to significantly penetrate into the sample. 
For strong disorder, 7 ::^ A, states in the surface layer 
become strongly localized, but new weakly disordered 
metallic states appear which reside primarily in the first 
inward layer of the sample. Remarkably, these states dis- 
play a clear Dirac dispersion, suggesting that they can be 
interpreted as metallic states at the interface between a 
TI and a strong Anderson insulator (AI). 

Our calculations demonstrate that already moderate 
disorder on TI surfaces requires to consider physics be- 
yond effective 2d Dirac models. On the application side, 
we propose that strong disorder can be used for nanos- 
tructuring of TI devices: Disordered layers in a TI crystal 
act as atomic-scale barriers which cut the TI. 

Model and Methods. As a minimal model for a 3d TI we 
take four bands in a cubic lattice of N sites, which may 
also be considered as an effective model for strained 3d 
HgTe and for insulators of the Bi2Se3 family [6|, Il5[ll6|. 
The Hamiltonian can conveniently be expressed in terms 
of the identity, F^ and four Dirac matrices F^, r*^^'^'^'"^^ = 
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Figure 1: (Color online) Spectral function A{k^uj)^ Eq. (|2]), along a path in the 3d Brillouin zone, obtained by ED for the TI 
model ([T]) on a 16 X 16 X 33 lattice with open boundaries. Histogram binning has been employed on the uo axis, with bin width 
t/50. The strength of the surface disorder 7 increases from a) to d). In the clean case, 7 = 0, bulk states exist for |c<;| > t, 
while surface Dirac cones appear near surface momenta (0,0), as is nicely visible e.g. near k = V 18]. (The small gap at V 
is a finite-size effect, while the "shadow" bands originate from the open boundary conditions.) The surface Dirac cones are 
destroyed for moderate disorder, 7 = 5t, but re-appear for large disorder, 7 = 20t. 
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matrices refer to orbital (spin) space: 
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^n is a four-component spinor at site n, t is the hop- 
ping amplitude and m a parameter which tunes the band 
structure. The model describes a weak TI with two Dirac 
cones per surface for \m\ < t, a strong TI with a single 
Dirac cone per surface for t < |m| < 3t, and a conven- 
tional band insulator for \m\ > 3t. Our numerical results 
below are obtained for m = 2t^ where the bulk band gap 
is A = 2t. 

We implement quenched disorder via random on-site 
energies e^ drawn from a uniform box distribution with 
width 7. For most of the paper, we will be concerned 
with surface disorder, i.e., Cn G [—7/2,7/2] for sites n on 
each surface of the sample, and e^ = otherwise. 

Using state-of-the-art exact diagonalization (ED) rou- 
tines [17| we calculate {|m)}, the two- fold Kramers de- 
generate eigenstates oil-L ^. Those can be visualized in 
both momentum and position space, via the momentum- 
resolved spectral function 

4 AN 

A(fc,a;) = ^^ Km|V(fc,i))|'5(w-w™) (2) 

j — \ m=l 

and the local density of states (LDOS) 

4 AN 

Pn{^) = '^Y1 l(^lV^(^n, j))l^^(^ - ^m) • (3) 

j=l m— 1 



e„> '|V'(fc,i')) = exp(i(fc|eA ))Sjj' is a Bloch 



state and {cn^^ ej |'^(rn',/)) = ^nn'^jj' a Wannier 

state, where \en ) and |e • ) denote the canonical basis 
vectors of position and bandindex space, respectively. 

Momentum-resolved spectral function. In Fig. [U we 
start by presenting results for A(k^uj) along a path fol- 
lowing the high-symmetry directions of the hulk Brillouin 
zone. In the clean case, Fig.[TJi, bulk states exist at ener- 
gies \uj\ > t, and sub-gap surface states are clearly visible: 
For the model ([T]) with m = 2t the Dirac cone is located 
at the surface momentum Tg = (0, 0) and consequently 
shows up in A{k.,uj) at 3d momenta k which equal T.^ af- 
ter projection onto the respective surface, e.g., at T [18|. 

While for weak surface disorder. Fig. [TJd, the spec- 
tral function has changed rather little (apart from minor 
broadening effects), the picture is drastically altered for 
moderate disorder in Fig.[lh, where for 7 = 5t the typical 
en is somewhat larger than the bulk gap. As expected, 
the bulk states are well defined, but any sub-gap states 
have lost their momentum-space structure. As will be- 
come clear below, the surface states have acquired diffu- 
sive character. Finally, Fig.[l]i is astonishing, as it shows 
that sub-gap states with a well-defined Dirac dispersion 
re-appear for strong surface disorder, 7 = 20t. 

Local density of states. More information can be ex- 
tracted from the LDOS Pn(<^) [211 • Fig. [2] displays the 
LDOS for the two outer layers (surface and first inward 
layer) of a finite system at cj = 0, i.e., in the middle of 
the bulk gap [22|. pn{oj) can further be used to construct 
the LDOS probability distribution f[pn{oo)] [23], i.e., a 
histogram of LDOS values at a fixed energy cj, for spe- 
cific parts of the sample. Fig.[3]shows such layer-resolved 
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Figure 2: (Color online) Spatially resolved LDOS Q in the middle of the bulk gap, pn{oj=^0), shown for the surface (top row) 
and the first inward layer (bottom row) of a 16 x 16 x 33 lattice [23]. p is normalized to the mean DOS of an extended bulk state, 
Pme = 1/(4A/'), such that the weight of surface states scales as N-^^^. Each panel shows the 6 faces according to the scheme 
in the lower right corner; the strength of the surface disorder 7 increases from left to right. While in the clean case, 7 = 0, 
surfaces states penetrate the bulk only slightly, this changes already for moderate disorder, 7 = 5t. For 7 = 20t tendencies of 
localization on the outer shell are visible, whereas weakly disordered states appear in the first inward shell. 



histograms for pn{uj = 0). 

In the clean case, Figs. [2^ and [3^, we have homoge- 
neous surface states, with the weight being located al- 
most exclusively in the outermost layer. These states are 
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Figure 3: (Color online) Distribution of the in-gap LDOS, 
Pn{(^ — 0)/y9me, for the surfacc and two next inward layers. 
Panels b,d,f) show different strengths 7 of surface disorder, 
and panels c,e) show bulk disorder, i.e., all tn G [—7/2,7/2] 
for comparison. While the clean case a) has sizeable low- 
energy LDOS only on the outer layer, weight is distributed to 
the inner layers with increasing disorder. For strong surface 
disorder, panel f), the broad distribution for the outer layer 
indicates localization, while disorder effects are much weaker 
on the high- intensity first inward layer. In contrast, for bulk 
disorder the broad LDOS distribution essentially coincides for 
all layers (surface states are destroyed) and indicates diffusive 
metallic behavior in c) and localization in f). 



only weakly perturbed for 7 = t, Figs.[2]3 and[3]3, consis- 
tent with topological protection (although the weight on 
the first inward layer is seen to increase). 

Moderate surface disorder, 7 = 5t, significantly re- 
distributes part of the weight of the in-gap states from 
the surface to the first inward layer (and even to the 
second inward layer). Figs. ^ and [3]i. Moreover, the 
states become inhomogeneous, and the LDOS distribu- 
tion significantly broadens. These are clear signatures 
of diffusive metallic behavior - this appears consistent 
with the fact that a 3d TI under influence of moderate 
hulk disorder becomes a diffusive metal [19|, |20|. (The 
latter is nicely visible in Fig. [Sh, which shows disordered 
low-energy states in all layers for a bulk 7 = bt.) We con- 
clude that, for surface disorder, low-energy surface states 
do exist for 7 = 5t, but their momentum-space Dirac 
structure is lost, consistent with Fig. [it. 

We now turn to strong surface disorder, 7 = 20t. Here, 
the main weight of the low-energy states has been pushed 
into the sample (!) to the first inward layer. Fig. [3f. In 
contrast, the LDOS in the surface layer is small with 
a very broad distribution, corresponding to the "frag- 
mented" real-space structure in Fig. [2]i. The surface 
layer thus shows a clear tendency towards strong An- 
derson localization. (This conclusion is again supported 
by the corresponding result for bulk disorder. Fig. [3^, 
signalling bulk AI behavior [24].) At the same time, the 
first inward layer displays a much narrower LDOS dis- 
tribution, i.e., is only weakly affected by disorder (recall 
that the Cn are zero here). Together with the data in 
Fig. [T]i, this shows that metallic low-energy states have 
emerged in the first inward layer and exhibit a well- 
defined Dirac structure in momentum space, much like 
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Figure 4: (Color online) Spatially resolved LDOS at energies cj = and t for al6xl6x33 sample with both surface and 
midplane disorder, for disorder strengths 7 = t (top) and 7 = 20t (bottom). The individual panels show the LDOS slice- by-slice 
for the entire sample. While weak disorder leaves the physics of the clean sample untouched, strong disorder leads to localization 
on the surfaces and in the midplane, and new weakly disordered sub-gap states appear both in layers 2 and 32 and in layers 16 
and 18 adjacent to the midplane: The disordered midplane has cut the TI into two. 



the surface states of the clean sample. The situation thus 
can be interpreted as an interface between a TI and an 
AI (which here consists of the surface layer only), with 
TI interface states forming in the outermost layer of the 
clean TI. 

In the limit of large 7, the behavior can be rational- 
ized by noting that the surface sites with le^l ^ t effec- 
tively decouple from the bulk, their only effect being to 
introduce an effective disorder in the first inward layer of 
strength t^/7, i.e., a smaller TI with weak surface disor- 
der obtains. 

Relation to surface Dirac theories. Given the drastic 
effects of disorder on TI surface states, a comment on 
disorder effects in 2d Dirac theories is in order. Recent 
work Ill4l3( has shown that even strong disorder cannot 
localize electrons of a single Dirac cone, which has lead 
to the notion of a supermetal and self-organized criti- 
cality in the presence of interactions. However, these 
theories are constructed in the continuum such that the 
cutoff energy of the Dirac spectrum is the largest energy 
scale. Therefore lattice effects are not captured, and it is 
precisely lattice effects which are responsible for strong 
localization in the limit of strong on-site disorder. (For 
TIs, such lattice effects include the coupling to sites away 
from the surface layer.) 

We conclude that the applicability of 2d Dirac theories 
to TI surface states is limited to energy scales small com- 
pared to the bulk gap A, and thus our results are not in 
contradiction to those of Refs.lll|,ll2|. Parenthetically, we 
note that attempts to study 2d lattice models with a sin- 
gle Dirac cone are usually hampered by the fermion dou- 
bling problem which can only be circumvented in special 
cases with broken time-reversal symmetry [25|. There- 
fore, treating the full 3d model as done here is indispen- 
sible. 

Midplane disorder: How to cut a TI. We can use the 
insights gained above for an "experiment": If a single 



strongly disordered layer, simulating an AI, produces new 
Dirac interfaces states in the TI, this can be used to effec- 
tively cut a TI. To this end, we introduce on-site disorder 
selectively in one central layer (or midplane) of the TI, 
Fig. m If this disorder is strong, 7 = 20t, the midplane 
shows signatures of Anderson localization, and new low- 
energy states appear to the left and right of the midplane. 
From Fig. [H one can deduce that these layers show the 
same characteristics as the layers adjacent to the disor- 
dered surfaces. 

Our simulated sample thus represents a layered struc- 
ture AI-TI-AI-TI-AI (the two outer AI on the sample 
surface), with interface states bounding both pieces of 
TI. (We have checked that the properties remain robust 
if the number of disordered layers is increased.) This 
suggests efficient nanostructuring via disorder: A large 
TI crystal can be "cut" into multiple TIs, preserving the 
single-crystalline structure, by intentionally disordering a 
few atomic layers, e.g., via irradiation, which then form 
barriers with new interface states. 

Conclusions. By large-scale numerics we have demon- 
strated that surface disorder can strongly disturb surface 
states of 3d TIs, such that their description in terms of 
a 2d Dirac theory breaks down: Topological protection 
is only relevant for weak disorder, while disorder with a 
strength comparable to or larger than the bulk gap first 
leads to diffusive metallic behavior and then to Ander- 
son localization at the surface. In the latter case, the 
emerging interface between a TI and an AI comes with 
new low-energy interface states, with properties similar 
to the surface states of a clean TI interfacing vacuum - 
the physics of these interface states clearly deserves fur- 
ther studies. 

STM experiments have already detected strong pertur- 
bations of surface states by isolated impurities [26| , and 
experiments on more disordered surfaces can be used to 
confirm our predictions. 
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